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Abstract. We present a self-dual bending theory that may enable a better 
understanding of highly nonlinear global behavior observed in biological vesicles. 
Adopting this topological approach for spherical vesicles of revolution allows us to 
describe them as frustrated sine-Gordon kinks. Finally, to illustrate an application 
of our results, we consider a spherical vesicle globally distorted by two polar latex 
beads. 
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1. Introduction 

Our primary motivation is to contribute to the understanding of highly nonlinear 
global behavior observed in vesicles. Such cases include closed vesicles with low 
genus exhibiting spontaneous conformal transformation [Q [21131; spherical vesicles 
globally distorted by a single latex bead 0], and spherical vesicles with several 
arms |5J- Whereas the difference in scale between the vesicle's membrane thickness 
and its overall size allows for the vesicle to be described as an embedded surface, 
the understanding of vesicle morphology is currently founded on the concept of 
bending elasticity [HI 13 El El EH ^2- To investigate vesicles, biophysicists widely 
invoke a harmonic bending Hamiltonian, historically introduced by Hclfrich and 
perform polynomial approximations to describe matter interactions. However, within 
the context of exploring nonlinear global behaviors such an approach may appear 
inappropriate to a geometric topologist. In this paper an attempt to construct a 
more suitable covariant field theory for the bending behavior of vesicles, outlined in a 
previous work ^2]; is illustrated through spherical vesicles of revolution. 

Customarily, the Monge representation (surface equations) of the vesicle shape 
leads to a local characterization of the surface by its two principal curvatures 
|13lll4lfT5| . and thereby to an expression for the bending energy as an expansion 
in the principal curvatures up to a given order and with respect to some desired 
symmetries: (i) the even symmetric expansion up to second order corresponds to the 
curvature energy density introduced by Canham [§] which has one phenomcnological 
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parameter; (ii) the most general symmetric expansion up to second order is the 
well known Hclfrich curvature energy density which contains three phenomcnological 
parameters |Sl[71|H1^2i (pfy the antisymmetric expansion up to second order fits the 
deviatoric bending contribution density suggested by Fischer |16U17| which has two 
phenomcnological parameters; (iv) expansion of higher order can be envisaged likewise 
|18| . By assumption, the general symmetric expansions in the principal curvatures can 
be formulated as expansions in the mean curvature (mean of the principal curvatures) 
and in the Gaussian curvature (product of the principal curvatures): the harmonic 
Helfrich curvature energy density is quadratic in the mean curvature and linear in the 
Gaussian curvature [H1I7II51ITT] . Commonly, as long as the vesicle remains in the same 
topological class, the Gaussian curvature term is dropped out jBHUIH] by invoking the 
Gauss-Bonnet theorem (13114115) which claims that the total integral of the Gaussian 
curvature over the vesicle surface depends only on the topology of the shape — the 
total Gaussian curvature measures the genus of the shape. From the perspective of 
a geometric topologist, ignoring such a global scale property may be undesirable and 
may essentially favor the notion of bending elasticity on the local scale. Although this 
local approach allows a large variety of phenomena to be understood OEMZUHlElEni , 
it may certainly preclude a deep understanding of some global bending phenomena as 
it fails to reveal the topological classification of vesicles. 

In contrast, the fundamental theorem of surface theory |13H14lH51l21| leads to 
a description of the surface by a prescribed metric tensor and a prescribed shape 
tensor related to each other by integrability conditions. We should therefore inquire 
how bending elasticity may be formulated within a covariant field theory for vesicles, 
and ultimately determine which principle may dictate bending elasticity. With this 
in mind, the observation of closed vesicles with low genus exhibiting spontaneous 
conformal transformation 02113] strongly suggests imposing a covariant functional 
which reveals the underlying topology and which is globally invariant under smooth 
conformal transformations as a bending Hamiltonian. In general, symmetries and 
conservation laws are connected by the Noether theorem |22M23ll2"4*) . Nevertheless, 
in this context, the conserved entity emerges from topology while the pertinent 
transformations are restricted by topology: two vesicles of different genus cannot 
be continuously deformable into one another — they are topologically distinct. As a 
matter of fact, topology gives rise to new physics [21121] with some of the following 
relevant features: (i) metastable configurations (mostly solitons) fall into distinct 
classes, of which the trivial configuration (vacuum) belongs only to one class; (ii) there 
exists no superposition principle, i.e. resultant configurations form complicated ones; 
(iii) frustration phenomena may arise from topological obstruction leading to both 
global and localized effects |2fill27ll2"£ll2"5] . 

Nccthcr-like calculational machineries exist to study topological systems: vesicles 
were shown ^2] to be subject to a technique known as the Bogomol'nyi decomposition 
which successfully treats various topological models in fields ranging from condensed 
matter physics to high energy physics [^OOIEIEIIES] ■ Concisely, the Bogomol'nyi 
technique applies to the total integral of the contracted self-product of the shape 
tensor, known to be globally invariant under conformal transformations (341 135) . 
revealing the topological nature of bending phenomena by both identifying the non- 
trivial metastable bending configurations and splitting bending configurations into 
topological classes according to their genus/end. The non-trivial metastable bending 
configurations are the round sphere and the minimal surfaces up to a conformal 
transformation of the ambient space (e.g. catenoid and Lawson surfaces j2E]) in 
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agreement with observations EI EI El E3| - Furthermore, it is worth noticing that 
any deformation of the non-trivial metastable bending configurations spontaneously 
matches, for vesicles of spherical topology, the deviatoric bending contribution 
described by Fischer |16U17j and, for vesicles of non-spherical topology, the mean 
curvature bending contribution (up to a conformal transformation of the ambient 
space) described by Helfrich IHIITIIHIITT] . For completeness, we note that the involved 
total integral corresponds to the covariant form of the bending energy proposed by 
Canham (JJj. Clearly, at least to study global bending phenomena, the topological 
approach is very appealing. Nevertheless, the mathematics involved here may 
discourage some readers. To overcome this, we show how this unorthodox approach 
leads to a description of spherical vesicles of revolution in terms of 'frustrated' sine- 
Gordon kinks. 



2. Self-dual approach 

Before focusing on spherical vesicles of revolution, we first succinctly expose how 
topology in vesicles is revealed by the Bogomol'nyi decomposition |24ll3()ll3*lll32ll33j . 
Following the fundamental theorem of surface theory, we represent the vesicle shape 
by a pair of tensors coupled to each other by integrability conditions, namely a 
prescribed first fundamental form (the metric tensor) coupled to a prescribed 
second fundamental form (the shape tensor) bij with respect to the equation of Gauss 
and the equations of Codazzi and Peterson |13U14l lT5| . For the bending Hamiltonian, 
we consider the covariant functional 

n b [s} = ^kfdsb l3 b^, (i) 

Js 

which depends on the vesicle shape S through the prescribed pair (g%j,bij) and on 
the phenomenological parameter k describing the bending rigidity. The summation 
convention has been adopted while customary notations have been used: the integral 
runs over the surface manifold S with surface element dS = da; 2 ^/[gf, where \g\ 
represents the determinant det(gij) and x the set of arbitrary intrinsic coordinates. 
Observe that the suggested covariant functional (QJ corresponds to the covariant form 
of the bending energy invoked by Canham [§] since the scalar bijb 1 * is equal to the 
trace of the square of the mixed shape tensor 6/ , and thus to the sum of its squared 
eigenvalues, namely the sum of the squared principal curvatures. Next, we define the 
dual tensor a^- of an arbitrary tensor by 

cUj = e ik eji a kl (2) 

where e mn denotes the totally antisymmetric tensor associated with the surface S. It 
is easily checked that the dual transformation <Zjj — ► Sjj is an involution, i.e. 

o%j — Q>ij • (3) 

While there exists no such relation between the shape tensor bij and its dual bij, the 
metric tensor gij satisfies the property of self-duality; that is 

9ij=9ij- (4) 

Nevertheless the bending Hamiltonian Q remains locally invariant under the dual 
transformation as we have, after straightforward index rearrangements, 

/>,,./<'•' = hjb ij . (5) 
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Most significantly, the extrinsic curvature G, which yields 

G=\b ij b ij , (6) 

remains clearly unchanged under the dual transformation according to the involutivc 
relation The interesting feature is that the extrinsic curvature G coincides with 
the intrinsic curvature (the Gaussian curvature) K when the embedding ambient space 
is flat: its total integral then measures the topology of the vesicle shape. 
More precisely, the Gauss-Bonnet theorem |13II38| claims that 



f dSG = -4tt (g + e-1) 



(7) 



where <S gje is a surface manifold embedded in a flat ambient space and topologically 
equivalent to a closed surface manifold of genus g less e points (ends). These 
considerations lead to the construction of a self-dual theory by introducing the self- 
dual/anti-self-dual tensors 



R± 1 



bij dz bij 



(8) 



which manifestly verify the property of sclf-duality/anti-sclf-duality according to the 
involutivc relation 

B± ij =±B ± ij . (9) 

Contracting the self-product of the self-dual/anti-self-dual tensors (JSJ and using the 
equalities (JSJ and © readily show that the bending Hamiltonian density Jtf, extracted 
from the bending Hamiltonian decomposes as 

J% = \k B ± ij B ±ij + kG. (10) 

Integrating this decomposition (|1C1|> over the surface manifold S and recognizing the 
total curvature Q leads to rewriting the bending Hamiltonian in each topological 
class specified by the pair (g, e) in the form 



U b [S^} = \k\ dSB ± lJ B ±lJ ±4^ (g + e-1). (11) 



is,, 

Since the tensors B^ij yield the precious inequalities 

B ± ij B ±ij > 0, (12) 

the decompositions (|llf> saturate when the tensors B ± ij vanish. Henceforth, since the 
shape tensor is self-dual (B~ij =0) only for the round sphere S 2 , the decomposition 
with sign (— ) in l|ll|) is relevant only for the vesicle surfaces So topologically equivalent 
to the round sphere S 2 ; we get 

H b [S ] = \k\ dS B- i5 B- ij + Ank. (13) 

Similarly, as the shape tensor is anti-self-dual {B + ij = 0) only when the surface 
manifold is a minimal surface, the decomposition with sign (+ ) in (|llfl is pertinent 
only for the vesicle surfaces <S g . e topologically equivalent to a minimal surface of genus 
g with e ends. Since within flat space there is no closed minimal surface (e = 0) whereas 
minimal surfaces of genus g (g^O) with e ends (e^2) do exist [HE], we write 



H b [5 g , e ] = ifc/dS ir.jr" + 4nk (g + e-1 



-1) (14) 
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which is valid only when g and e ^ 2. Notice that the minimal surface of genus 
zero with two ends is the catenoid. Tedious considerations allow one to expand the 
previous decomposition (|14|) to closed surfaces of arbitrary genus (g 1) ^21- Besides 
bringing out the inherent topology, the Bogomol'nyi decompositions (|13f) and 1 1141 
offer a new perspective on bending phenomena: when the shape tensor obeys the local 
property of self-duality/anti-self-duality then the bending Hamiltonian is saturated. 
Conversely, when the shape tensor violates the local property of self-duality /anti-self- 
duality then the bending Hamiltonian is frustrated, i.e. an extra energy contribution 
is spontaneously created [1211261127112811291 that tends to vanish globally in the system. 
Before following this matter any further, we turn our attention to the spherical vesicles 
of revolution which are of pertinent interest. 

3. Spherical vesicles of revolution 

Now, we concentrate on vesicle shapes that are both smoothly transformable to the 
round sphere § 2 and can be generated by rotating a two-dimensional curve (a profile) V 
about an axis, i.e. on spherical vesicles of revolution. For our purpose, the azimuthal 
isometric coordinates (u, ip) appear to be a suitable choice of coordinates as follows. 
First, adopting isometric coordinates ('isothermal' coordinates) allows us to write the 
metric tensor g^ in the form |13ll25ll59"] 

9uu = 9 VV = e 2ff(u) (15) 

where the local Weyl gauge field a @U| depends only on the isometric coordinate u 
because of the azimuthal symmetry. Moreover, for any surface of revolution, analytical 
manipulations show that the shape tensor bij in azimuthal isometric coordinates (u, ip) 
yields 

b uu = ~e a{u) d u n(u) and b lpip = ~e' j{u) sinfi(it), (16) 

where Q corresponds to the polar angle of the unit normal vector of the surface. 
However, in order for the two tensors g^ and bij, cast respectively in the forms (11511 
and l|16|) . to be the first and second fundamental forms for a surface 5 in I 3 , they 
must satisfy the Gauss equation and the Codazzi and Peterson equations |13U14llT5| : 
in this case, after easy computations, this set of integrability conditions reduces to the 
equation 

d u a(u) = cos Q(u). (17) 

Conversely, simple calculations show that the Gauss- Wcingartcn equations |13lll4llT5] 
associated to the first form l|15|l and the second form (|16|) under the integrability 
condition l|17|) give a unique surface of revolution, except for both its position and 
scale in space. To summarize, any pair of tensors {gij, bij) that verifies the formulae 
11511 . (|16fl and l|17fl in azimuthal isometric coordinates (u, ip) is equivalent to a surface 
of revolution S in M 3 . Substituting H15fl and (|16f) into © leads to 

B ± uu = ±B ± VV = =fe e" (u) [d u n(u) ± sin n(u)] . (18) 

Henceforth, by applying (|13|) . the bending Hamiltonian Bogomol'nyi decomposition 
for the spherical surfaces of revolution So reads 

H b [S ] =Tik[du [d u n(u) -smft(u)] 2 +4tt/c, (19) 
Jv 
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where the integral runs along the profile V of the spherical surface So- The bending 
Hamiltonian i|19|) is saturated when the polar angle f2 satisfies the equation 

d a Q(u) =sinQ(u), (20) 

whose centered solution is the sine- Gordon kink 

O(u) = 7r — arccos(tanh(M)). (21) 

Resolving the Gauss- Weingarten equations |13U14llT3] with respect to (II 51) . (11(51) . (II 711 

and l|21f) uniquely gives (except for its position and scale in space) the surface, in 
cylindrical parametrization, 

r(u) = sech(u), z{u) = — tanh(w), (22) 

which is the round sphere § 2 in azimuthal isometric coordinates (it, <p), u running from 
— oo to +oo. Hence, in accordance with the existence theorem previously invoked, the 
spherical vesicle of revolution which saturates the bending Hamiltonian is the round 
sphere S 2 . However the key result is contained in the fact that the mctastablc spherical 
vesicle of revolution can be described by a saturated sine-Gordon kink: any spherical 
vesicle of revolution may be thereby envisioned as a frustrated (i.e. unsaturated) 
sine-Gordon kink, and highly nonlinear global physics may henceforth emerge from 
the model |26I27I28I29 |. Therefore, without loss of generality, the initial bidimensional 
system has been reduced to a well known unidimensional nonlinear system. Next we 
illustrate our approach through an illuminating example. 



4. Latex bead on a spherical vesicle 



Grafting a latex bead on a spherical vesicle induces a global distortion characterized by 
two polar concave regions: a pinched concave region around the bead faces a smooth 
one as if a second bead were diametrically opposite @]. Accordingly, we assume a 
bare spherical vesicle of revolution with two identical round latex beads which are 
respectively attached at the north and south poles. By bare, we mean that only the 
bending Hamiltonian is considered. To begin with, substituting (|15fl and into 
JQ) gives the Hamiltonian of our system in the canonical form 



H b [S] = irk / du d u Ct(u) + sin 2 Vt(u) 



(23) 




Figure 1. North pole sketch of a spherical vesicle of revolution distorted by two 
identical round latex beads grafted at its poles as deduced by the self-dual bending 
theory. The dashed circle profiles the bead, the bold arc the polar cap imposed 
by it. The arrows indicate the nomenclature: the latex bead dimensionless 
radius p, the encapsulated dimensionless radius p, and the encapsulation angle 
a. Dimensionless latex bead radius p = 0.15; relative encapsulated radius p/p= ^. 




Figure 2. North hemispheric profiles of a spherical vesicle of revolution distorted 

by two identical round latex beads grafted at its poles with respect to different 

encapsulated radii p as deduced by the self-dual bending theory. The bold arcs 

correspond to the profiles of the polar cap imposed by the bead. Dimensionless 

latex bead radius p = 0.15; relative encapsulated radius p/p from outside to inside: 
n 1 3 1 5 3 r i 

u > 4 > 8 ' 2 ' 8 ' 4 > 8 > 



The Euler-Lagrange equation derived from (|23|) (or from is the sine-Gordon 

equation 

d uu n(u) = sinf2(u) cosf2(u). (24) 
The solution of J23J) extending (I2TT) . 

Q(u) — 7r - arccos(sn(^= | m)) m G (0, 1] , (25) 
corresponds to the surface of revolution 



r{u) = t+7 
z ( u ) = 7T77 



dn( _^| m) + ^ cn( _^| m) 
msn(-^|m)+Dn(J=|m) 



( 1 — m ) u 



(26) 



which smoothly joinsf the concavely bounded spherical polar caps to the vesicle surface 
when the parameter m satisfies 

2 



(27) 



_p(l + p 2 ) + 2p 2 _ 

p and p being respectively the latex bead and encapsulated dimensionless radii — 
see Figure Q] and Figure [21 The isometric coordinate u runs from — u to +u with 
u = \/m[2K(m) — F(cosd | m)], where the encapsulation angle a yields sin a = p/p. 
The reader is encouraged to check that, when the encapsulated dimensionless radius p 
vanishes, the parameter m effectively tends towards 1 and thus the surface (|2^|) tends 
to the round sphere S 2 if^ as expected. Now the extra energy AE spontaneously 

f The junction conditions are the ones implicitly assumed in Ref. Q]: at the junction points, the 
normal lines of the vesicle surface and of the latex bead surface are imposed to coincide — see FigurclTl 
Such an assumption is reasonable at the vesicle scale, nevertheless at the membrane scale an adequate 
treatment may be needed since the curvature experiences a discontinuity at the junction points — as 
far as we know there is no such treatment in the literature. 
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Figure 3. The relative extra bending energy AE/4nk versus the relative 
encapsulated radius p/p for different radii p of the grafted latex beads: the radius 
p increases in the direction of the arrow. Dimcnsionlcss latex bead radius from 
right to left: 0.01, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. 



generated by the polar latex beads can be computed. We have 

-l 

<p<p 73 



B uu B Vlf> — -k j— (28) 



from which we obtain 
4-rrk " 



AE 



2 [E(m) - K(m)] - [E(cos a \ m) - ^ F(cos a | m)] • 



.(29) 



Standard notations have been adopted! jlllH2lHHj . Observing that, according to 
l|27ll . the parameter m decreases strictly from 1 to with respect to the encapsulated 
radius p, allows an effortless description: when the encapsulated radius p vanishes, 
the underlying soliton (|25fl reaches its saturated configuration (|21|) ; when the 
encapsulated radius p increases, the underlying soliton (|25(l becomes frustrated and 
the bending energy is accordingly altered — see Figure Finally, our naive model 
system demonstrates at least one weakness of the traditional (local) approach for the 
description of global bending phenomena, as follows: since its mean curvature H, or 
the half trace of the shape tensor, is constant, 

H=\^i (30) 

the surface of revolution (|26|) causes the harmonic Hclfrich bending Hamiltonian 
density either to vanish if a phcnomcnological spontaneous curvature is introduced or, 
if not, to be proportional to the square of the mean curvature PI7I I81HT] . irrespective 
of the value of the encapsulated radius p. Consequently the traditional approach 
is partially insensitive to the encapsulation mechanism experienced by our system 

| The functions sn, cn and dn are the Jacobian elliptic functions, F and K the incomplete and 
complete elliptic integrals of the first kind, E the incomplete and complete (depending on the 
argument) elliptic integrals of the second kind, and the function Dn denotes the Jacobian Epsilon 
function defined as Dn(ii | m) = f^At dn 2 (t \ m) 1411 . 
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in the sense that an observer somewhere on the surface outside the grafted regions 
cannot determine whether or not latex beads are grafted at the poles just from 
the harmonic Hclfrich bending Hamiltonian density, whereas the same observer can 
obtain a quantitative answer by computing the extra bending Hamiltonian density. 
In the former case, what occurs globally is not known because in every scenario the 
harmonic Hclfrich bending Hamiltonian density takes the same form; in the latter, 
the extra bending Hamiltonian density vanishes when there are no latex beads but 
it measures the deformation otherwise. Nevertheless it may be objected that, if 
no phenomenological spontaneous curvature is introduced, computing the harmonic 
Hclfrich bending Hamiltonian density allows for an answer which is quantitatively 
similar, provided that the expression is known for the particular state, the one without 
latex beads for example. Such exact information, required for the former approach 
but superfluous for the topological approach, may not be so easily identified in a more 
realistic case. Second, the topological approach provides actually a finer response as 
the deformation measurement is really enabled by the anti-sclf-dual tensor which can 
be conceived as a bending deformation tensor. 

5. Conclusion 

In conclusion, the results of our investigations not only provide a novel viewpoint 
regarding bending phenomena, but also demonstrate the possibility of describing 
spherical vesicles of revolution by frustrated sine-Gordon kinks as a method to 
investigate the essence of highly nonlinear global behavior observed in vesicles. The 
geometry of vesicles has been studied quite extensively in recent years |44I45I46I47| but 
it was not considered in conjunction with the global topological aspects. Of course, the 
present model is too naive to capture the complexity of biological vesicles, in particular 
it should be complemented by the inclusion of material fields and constraints H3EIIE] . 
Still, as far as bending phenomena are concerned, the self-dual bending theory should 
provide a powerful theoretical framework to study global behavior observed in vesicles 
or, in general, other deformablc systems. 
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